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Abstract. This paper defines a class of left invariant variational problems on 
a Lie group G whose Lie algebra g admits Cartan decomposition g = p + t with 
the usual Lie algebraic conditions 

[p,p] CI , [p,l] cp, [i,t]ct 

The Maximum Principle of optimal control leads to the Hamiltonians H on 
g that admit spectral parameter representations with important contributions 
to the theory of integrable Hamiltonian systems. Particular cases provide nat- 
ural explanations for the classical results of Fock and Moser linking Kepler's 
problem to the geodesies on spaces of constant curvature and J.Moser's work 
on integrability based on isospectral methods in which C. Newmann's mechan- 
ical problem on the sphere and C. L. Jacobi's geodesic problem on an ellipsoid 
play the central role. The paper also shows the relevance of this class of 
Hamiltonians to the elastic curves on spaces of constant curvature. 



1. Introduction 



A Lie group G with an involutive automorphism a admits several natural vari- 
ational problems whose solutions provide new insights into the theory of integrable 
Hamiltonian systems and to the geometry of the associated homogeneous spaces. 
An involutive automorphism ct on a Lie group G induces a splitting g = p © { of 
the Lie algebra g of G with { equal to the Lie algebra of the group K of fixed points 
under a. When G is semisimple then p is the orthogonal complement to t relative 
to the Killing form and p and 6 satisfy the following Lie algebraic relations: 

(1-1) [p,P] = e, [p,«] =p- 

The first relation implies that any two points of G can be connected by a curve 
whose tangent takes values in the left invariant distribution ^{g) = {gU : U G p}. 
In the case that (G, K) is a Riemannian symmetric pair there is an AcIk invariant, 
positive definite quadratic form ( , ) on p that induces a natural optimal control 
problem on G : minimize the integral ^ /o^(^^(^) iC^(i) )dt among all curves 
g{t) G G that are the solutions of 

(1.2) ^=.g{t)U{t), U{t)ep, te[0,T] 
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with fixed boundary conditions g{0) ~ go ,g{T) = gi. Here go and gi are arbitrary 
but fixed points in G and the terminal time T > is also fixed. This problem, 
called the canonical sub-Riemannian problem on G, is well defined in the sense 
that for any pair of boundary points in G there is an optimal solution. It is well 
known that optimal solutions are of the form 

(1.3) g(t)=5oe*(^+Q)e-*0 

for some elements P E p, and Q G i ([17] ). The above implies that any element g in 
G can be represented as g = e^^+'^^e"'^ for some elements P S p, and Q € t. This 
sub-Riemannian problem is naturally related to the canonical Riemannian problem 
on the quotient space M = G/K in the sense that the Riemannian geodesies are 
the projections of the above curves with (5 = 0. 

In this paper we will be interested in another optimal control problem defined 
by an affine distribution 'D{g) = {g{A -\- U) : U £ t] with A a regular element in 
p under the assumption that the Killing form is definite on 6 ( which is true when 
if is a compact subgroup of G ) . This optimal control problem consists of finding 
the minimum of ^ {U (t) , U (t)) dt among aU solution curves g{t) e G of the 
affine control problem 



(1-4) ^{t)=g{t){A + U{t)),U{t)et, te[0,T] 

subject to the given boundary conditions 17(0) = go, g{T) = gi where the qua- 
dratic form ( , ) denotes a scalar multiple of the Killing form that is positive definite 
on I. 

This problem might be regarded as the canonical affine problem on symmetric 
pairs (G, K) for the following reasons. The reciprocal affine system with A £ t 
and L/ € p is isomorphic to (11.21) and bears little resemblance to the solutions of 
(|1.4p . Moreover, two affine systems defined by regular elements Ai and A2 are 
conjugate. 

The affine problem ()1.4|) will be referred to as the Affine-Killing problem or (Aff ) 
for brevity. It is first shown that (Aff) is well defined in the sense that for any pair 
of boundary conditions (50,51) there exist T > and a solution g{t) of (jl.4|) that 
satisfies 5(0) = 50, g[T) = 51 such that the the control U{t) that generates g{t) 
mmimizes the integral i {U{t), U{t))dt among all other controls whose solution 
curves satisfy the same boundary data. Then it is shown that optimal solutions 
are the projections of the integral curves of a certain Hamiltonian system on the 
cotangent bundle T*G oi G obtained through the use of the Maximum Principle 
of Optimal Control. 

To preserve the left invariant symmetries, T*G is realized as the product G x g* 
with 0* equal to the dual of the Lie algebra g of G and then g* is identified with 
via the Killing form so that ultimately T*G is realized as G x g. In this setting 
the Hamiltonian associated with (Aff) is of the form 

H ^^{L,,Lt) + {A,Lp) 

where Lp and denote the projections of i G onto the factors p and 6. 
Because K acts on p via the adjoint action, g as a vector space carries two 
Lie algebras: a Lie algebra of G and the Lie algebra g^ of the semidirect product 
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Gs ^ K >i p. The affine problem then admits an analogous formulation on Gs 
and the Maximum Principle leads to the Hamiltonian H that formally looks the 
same as the one obtained in the semisimple case. We refer to the semidirect version 
of (Aff) as the semidirect shadow problem. 

The essence of the paper lies in the integrability properties of the associated 
Hamiltonian flows H and Hs . It is shown that each flow admits a spectral repre- 
sentation 

(1-5) ^^[Ah,L,] 

where Mx — j{Lp — eA) and L\ — Lp — XLf + (A^ — s)A, with s = 1 in the 
semisimple case and s = in the semidirect case. Hence, the spectral invariants 
are constants of motion for the associated Hamiltonian flows. 

On spaces of constant curvature these results recover the integrability results 
associated with elastic curves and their mechanical counterparts ( [17;). Remark- 
ably, the spectral invariants above also recover the classical integrability results C. 
Newmann for mechanical systems with quadratic potential ('22') and the related 
results of C.G.J. Jacobi concerning the geodesies on an ellipsoid. The present 
formalism also clarifies the contributions of J. Moser ( [20]) on integrability of 
Hamiltonian systems based on isospectral methods. More significantly, this study 
reveals a large class of integrable Hamiltonian systems in which these classical ex- 
amples appear only as very particular cases. 

2. Notations and the Background material 

The basic setting is most naturally defined through the language of symmetric 
spaces. The essential ingredients are assembled below. 

An involutive automorphism cr on G is an analytic mapping G G, a ^ I that 
satisfies 

(2.1) ct(525i) = CT(.g2)cr(gi), for all gi, g2 in G. 

Then the tangent map cr* of a induces a splitting g = p J of the Lie algebra q 
of G with 

(2.2) p ={A e g : cr,(A) = -A} and k= {A <E g : a^{A) = A} 

The fact that ct* is a Lie algebra automorphism easily implies the following Lie 
algebraic relations 

(2.3) [P,P]C« , [p,f]Cp, [t.,t]Ct 

It follows that 1 is a Lie subalgebra of g, equal to the Lie algebra of the group 
K = {g (z G : (j{g) = 5} and that p is an AdK invariant vector subspace of g in the 
sense that Adh{p) C p for any h ^ K. An AdK invariant non-degenerate quadratic 
form on p will be called pseudo Riemannian. It is easy to show by differentiating 
that an AdK invariant quadratic form < , > on p is invariant, in the sense that 

(2.4) < A,[B,C]>=<[A,B],G > 

for any A, C in p and any B in ?. 

A pseudo Riemannian form that is positive definite will be called Riemannian. 
In the literature of symmetric spaces ([6]) the pair {G,K), with K a closed sub- 
group of G obtained by an involutive automorphism on G described above, is 
called a symmetric pair. If in addition this pair admits an AdK invariant positive 
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definite quadratic form ( , ) on p tlicn it is called a Riemannian symmetric pair. 
Riemannian symmetric pairs can be characterized as follows. 

Proposition 1. Let Adh,p denote the restriction of Adh to p. Then a symmetric 
pair (G, K) admits a Riemannian quadratic form { , ) on p if and only if 
{Adh,p : h € K} is a compact subgroup of Gl{p). 

In the text below we will make use of the Killing form {A, B)k = Tr{adAo adB) 
for A and B in g, where Tr{X) denotes the trace of a linear endomorphism X. The 
Killing form is invariant under any automorphism on g and in particular it is 
AdK and Ada invariant ([6]). The invariance relative to Ada implies 

(2.5) {A,[B,C])k = {[A,B],C])k 

for any matrices A,B,C in g. Spaces p and t are orthogonal relative to ( , )k 
because {A,B)k = {a^{A),a^,{B))k = {—A,B)k = — {A,B)k for any ^ in p and 
any B in t. 

In this paper (G, K) will be assumed a symmetric Riemannian pair with G 
semisimple and connected and K compact. Semisimplicity implies that the 
Killing form is non-degenerate, which then implies that its restriction to p is 
pseudo Riemannian. Semisimplicity also implies that the Cartan relations (j2.3p 
take on a stronger form 

(2.6) [P,p]-fi, [P,«]=P, 

The fact that K is a. compact subgroup of G implies that the Killing form is 
negative definite on 4 ([5], p. 56). In the sequel ( , ) will denote any scalar multiple 
of the Killing form which is positive definite on {. Under these conditions then \ \U\\ 
will denote the induced norm \\U\\ — (C/, U). 

An element A in p is said to be regular if {B e p : [A, i?] = 0} is an abelian 
algebra. It follows that A is regular if and only if the algebra A spanned by {B E 
p : [A, B] = 0} is a maximal abelian algebra in p that contains A (j5j ). 

With these notions at our disposal we return now to the affine problem defined 
above. It will be convenient to adopt the language of control theory and regard 
(jl.4p as a control system with U{t) playing the role of control. In order to meet the 
conditions of the Maximum Principle control functions are assumed bounded and 
measurable on compact intervals [0,r]. Solutions of ( 11.41 ) are called trajectories. 
A control U (t) is said to steer go to gi in T units of time if the corresponding tra- 
jectory g{t),tE[0,T] satisfies 5(0) = 50, 5(7") = 5i- A trajectory (7(i) generated 
by a control U (t) on an interval [0, T] is optimal relative to the boundary condi- 
tions (go,gi) if the integral ^ \\U{t)\\'^ dt is minimal among all other controls 
that steer go to gi in T units of time. Controls that result in optimal trajectories 
are called optimal. Thus every optimal control U{t) gives rise to a unique optimal 
trajectory because the initial point go is fixed. 

3. The existence of optimal solutions 

Proposition 2. // A is regular then (Aff) problem is well posed in the sense 
that for any pair of boundary conditions (.90,31) there exist a time T > and a 
solution g(t) on the interval [0,T] that is optimal relative to go and gi. 

The proof of this proposition requires several auxiliary facts from the optimal 
control theory and from the theory of symmetric spaces. We begin first with the 
facts from the theory of symmetric spaces ([5], [6]). 
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A Lie algebra g is said to be simple if it contains no ideals other that {0} and 
0. A Lie group G is said to be simple if its Lie algebra is simple. The first fact is 
given by 

Lemma 1. // (G, K) is a symmetric pair with G simple then Adx acts irreducibly 
on p . 

Proof. Let V denote an Adx invariant vector subspace of p. Denote by the 
orthogonal complement of V in p relative to the Killing form ( , )k- Since 

{V,Adh{V^))k = {Adh-i{V),V^)k, it follows that is also Adx invariant. 
Therefore, [I, V]CV and [I, V^] C V^, which in turn implies that [V, V^])k = 0. 
It follows that \V, V-^] = by semisimplicity of q. 

The above implies that V + \V, V] is an ideal in g that is orthogonal to . Since 
Q is simple V + [F, = q and therefore V-^ — 0. But then V — p. □ 

The other facts which are needed for the proof are assembled below in the forms 
of propositions. 

Proposition 3. Suppose that (G, K) is a symmetric Riemannain pair. There exist 
linear subspaces pi, p2, . . . , pm of p such that 

(1) p = pi® p2 ® ® pm and 

(2) pi, p2,..., pm oltc pairwise orthogonal relative to the Killing form. 

(3) Each pi is ad{t) invariant and contains no proper ad{t) invariant linear 
subspace. 

Proposition 4. Let Qi = pi +[ pi , pi ] i — I, . . . ,m in a semisimple Lie algebra 
Q. Then, (1) Each Qi is an ideal of g and also a simple Lie algebra. 
Moreover, 

(2) [ Qi, Bj] = and ( g,, Qj)k = 0, i =^ j, and {X e Q : [X,Y] ^ Q ,Y q) ^ o. 

Corollary 1. Let Ai denote the projection of a regular element A in p on pi. Then 
each Ai ^ 0,i — 1, . . . , m. 

Proof. If Ai were equal to then [A, pi] = by (2) in Proposition [3l This would 
imply that p; is abelian by regularity of A, which in turn would imply that Qi — pi. 
But that would contradict (3) in Proposition [3l □ 

We now turn attention to the pertinent ingredients from the accessibility theory 
of control systems. The Lie Saturate of a left invariant family of vector fields J- 
is the largest family of left invariant vector fields (in the sense of set inclusion ) 
that leaves the closure of the reachable sets of invariant (H^)- It is denoted by 
LS{T). 

Since left invariant vector fields are defined by their values at the identity, the 
Lie saturate admits a paraphrase in terms of the defining set P in g. For the affine 
system dH ), T = {A + B : B e i}. 

Definition 1. Let P C g . The reachable set of P denoted by A(T) is the set of 
terminal points g{T) G G corresponding to the absolutely continuous curves g(t) on 
intervals [0,T] such that g{0) = L and ^{t)g~^{t) G P for almost all t G [0,T]. 

Then LS'(P), the Lie Saturate of P can be described as the largest family in q 
such that 



(3.1) 



^(^(^^(P)) =d(yl(P), 
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where cl{X) denotes the topological closure of a set X. 

The following lemma is well known in control theory ([15)). 

Lemma 2. a. LS{T) = g is a necessary and sufficient condition that A{T) — G. 

If C denotes the convex cone spanned by "^^^i OLiAdhX^) i hi € 4, Ui ^ 0, i = 
1, 2, . . . , m, m e Z+ then 

(3.2) ClMQLS{V) 
where T {A + [/:[/ G fi}. 

With these results at our disposal let us turn to the proof of the proposition. 

Proof. Let Traj{go, gi) denote the set of solutions g{t) of (|1.4p that satisfy ^(0) — 
go, g{T) — gi for some T > 0. If Traj{go, gi) is not empty for any go and gi in 
G then (II. 4p is said to be controllable. An argument based on weak compactness 
of closed balls in Hilbert spaces shows that there is an optimal trajectory g{t) in 
Traj{go, gi) generated by a control U{t) in L^([0,T]) whenever Traj{go, gi) is not 
empty (Theorem 1 in ( 17 ). But then it can be shown that an optimal control 
in L^([0,T]) is absolutely continuous and hence belongs to L°°([0,T]). The above 
argument shows that controllability implies the existence of optimal trajectories. 

Now address the question of controllability. According to Lemma 2 it would 
sufRce to show that the convex cone C spanned by {^aiAdhi{A), hi € t, ai ^ 
} is equal to p. 

Let V denote the vector space spanned by {Adh{A) : h G i } and let Vi denote 
the vector space spanned by {Adh{Ai) : h G i } where Ai is the projection of A on 
pi, as in the Corollary above. Each Vi is a non-zero Adx invariant vector subspace 
of a simple Lie algebra g^. According to Lemma [1] Vi — pi,i = 1, . . . ,m and hence, 
y = p. It follows that C = {J2 ctiAdhi (A), hi ^ i, ^ } is an Adx invariant 
convex cone with a non empty interior in p. 

Then C = p if and only if the origin in p were contained in the interior of C. 
Let 5" = {a: G p : ||X|| = ||A||}. If were not in the the interior of C, then Cn S"" 
would be a convex cone in the sense of Eberlein ([5J, 1.15) that is invariant under 
AdK- But then the sole of this convex set would be a fixed point of Adx which is 
not possible since Adx acts irreducibly on each pj. □ 

3.1. Semidirect products and the shadow problem. Recall that if Kq is 
a Lie group which acts linearly on a finite dimensional vector space V then the 
semidirect product Gg = V Kq consists of points {v, h) in V x Kq with the group 
operation (wi, /ii)(w2, /i2) = ivi + hi{v2), hih2). The Lie algebra of consists of 
pairs {A, B)mVx to where io denotes the Lie algebra of Kq, with the Lie bracket 

Pi, Bi), (^2, Sals = (Bi(^2) -52(^i),[Bi,B2]). 

Every semidirect product V k Ko admits an involutive automorphism a{x, h) = 
(— x, h) for every (x, h) G V x Ko. It follows that K = {0} x Ko is the group of 
fixed points of a and that 

(3.3) p = y X {0} and « = {0} x Jo- 
lt is easy to check that Adh{x,Q) = {h{x),0) for every h G K and every x G 
V. Therefore, (G, K) is a symmetric Riemannian pair if and only if Ko is a 
compact subgroup of Gl{V). 

Every Lie group G that admits an involutive automorphism carries the semidi- 
rect product Gs — p X K because K acts linearly on the Cartan space p via the 
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transformation 

(3.4) h{A) = Adh{A), Aep, for each h€K. 
Therefore, the Lie bracket on gs is given by 

(3.5) [{Ai,Bi),{A2,B2)]s - iadBi{A2) ~ adB2iAi),[Bi,B2]). 

If {A, _B) in p x: J is identified with A + B in p + t then the semidirect Lie bracket 
[ , ]s can be redefined as 

(3.6) [{Ai + Bi), iA2 + B2)]s - [Bi,A2] - [B2, Ai] + [Si, B2], 
from which it folfows that 

(3.7) [p, p], - 0, [p, e], = [p, i], [i, i]s = [e, t]. 

Thus is the underlying vector space for both Lie algebras g and Qs , a fact which 
is important for the subsequent development. The passage from to g can 
be described by a continuous parameter s by deforming the Lie algebra g^ to g via 
the Lie bracket [ , ]s '■ 

(3.8) [(Ai + Si), {A2 + B2)]s = [Si, A2] - [S2, Ai] + [Si, S2] + s[Ai,A2] 

We now return briefly to the afhne problem Aff to note that the data which is 
required for its formulation on a semisimple Lie group G also permits a formulation 
on the semidirect product Gg. The semidirect version consists of minimizing the 
integral i J^{U{t), U{t)) dt over all solutions g{t) in Gg of 

^^{t)=g{t){A + U{t)),U{t)el, [0,r] 

that meet the boundary conditions 17(0) — go, g(T) — gi. This "shadow" problem 
will be referred to as (AfTg). The same arguments used in the semisimple case show 
that (AflFs) is also well defined in the sense of Proposition [2l 

4. Left invariant Hamiltonian systems and the Maximum Principle 

Consider now the necessary conditions of optimality provided by the Maximum 
Principle. The Maximum Principle states that each minimizer is the projection 
of an extremal curve in the cotangent bundle T*G and each extremal curve is an 
integral curve of a certain Hamiltonian vector field on T*G. To state all this in 
more detail requires additional notation and terminology. 

As already stated earlier, g* denotes the dual of g. The dual of a Lie algebra 
carries a Poisson structure inherited from the symplectic structure of T*G realized 
as the product G x g* via the left translations. Functions on g* are called left- 
invariant Hamiltonians. If / and h are left-invariant Hamiltonians then their Poisson 
bracket {/, h} is defined by {/, h}{l) — l{[df, dh]), for I e g. 

On semisimple Lie algebras g* can be identified with g via the quadratic form 
( , ) with (L, X) = 1{X) for all X G g. In this identification p* and I* are identified 
with p and t whenever g admits a Cartan decomposition g = p © The above then 
implies that I = Ip + h with Zp S p* and 1^ e V is identified with L = Lp + Lf 
where lp and If correspond to Lp G p and L( e t. 

To preserve the left invariant symmetries T*G will be trivialized by the left 
translations and considered as the product Gx g*. The advantage of the above 
choice of trivialization is that the Hamiltonian lift of a left invariant vector field 
becomes a linear function on g*. Recall that a Hamiltonian lift of a vector field 



8 



VELIMIR JURDJEVIC 



X on a manifold M is a function Hx on r*Af defined as Hx{£,) — £,{X{x)) for 
each E T*M ([ISJ. If ^(g) — gA is a left invariant vector field on a Lie group 
G, then Hx{g,l)) = 1{A), Igq*. 

Any left invariant function h generates a Hamiltonian vector field h on G x 
g* whose integral curves {g{t),l{t)) are the solutions of the following differential 
equations 

(4.1) I = git)dhm), I = ^ad*idhmiit)), 

where dh denotes the differential of h considered as a element of g under the 
isomorphism (g*)* < — > g, and where ad*{L) : g* ^ g* is defined by ad*{L){l) — 
I o ad{L). 

With these notations at our disposal we now apply the Maximum Principle to 
the affine problem. The affine problem defines "cost-extended system" in R x G: 

(4.2) § = Imm^ ^{t) = git)iA + Uit)), U{t) e t. 
The Hamiltonian lift of the cost-extended system is given by: 

(4.3) Hu{\l) = A^IIC/II^ +;(yl + C/),A e R, / e g* . 

The above is a function on r*(RxG) tivialized as (R x R) x (Gx g*) with coordi- 
nates {x, X,g, I). Each control function U{t) generates a time varying Hamiltonian 
i/[/(t)(A,Z) ; the integral curves ^(f) — {x{t), X{t), g{t),l{t)) of the associated Hamil- 
tonian vector field i?(7(t) are the solutions of 

OH 

It follows that A is constant for any solution ^(t) since — = 0. 

Proposition 5. The Maximum Principle. Assume that tj{t) is an opti- 
mal control that generates the trajectory g{t). Let x{t) denote its running cost 
Jo ^II^IP*^^- Then {x{t),g{t)) is the projection of an integral curve of ^{t) = 
{x{t), X,g{t),l{t)) of Hijfj.-^ that satisfies the following conditions: 

(4.5) A < 0. When A then, l{t) ^ 0. 

(4.6) Ho(t)(X~m) > Hu{X~l{t))),U (.l.a.e. in [0,r]. 

In the literature on optimal control it is customary to consider only the projec- 
tions {g{t),l{t)) of integral curves ^(t) of Hu{X,l) which are parametrized by a 
non-positive parameter A. Control functions U{t) on T*G are called extremal if 
they generate solutions of (|4.4I) that satisfy conditions (|4.5p and (|4.6p of the Max- 
imum Principle. Extremal curves that correspond to A = are called abnormal 
and those that correspond to A < are called normal. In the normal case A is 
reduced to —1 because of the homogeneity properties of Hu{X, I) with respect to A. 

The Maximum Principle can be restated in terms of the extremals by saying that 
each optimal trajectory is the projection of an extremal curve (normal or abnormal) . 
Thus the Maximum Principle identifies two distinct Hamiltonians associated with 
each optimal control problem depending on the value of A. 
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Return now to the Hamiltonians of the affine problem. After the identifications 
oi I = Ip + It with L = Lp + Lt the Hamiltonian Hu{X,l) is identified with 
H{X,L)^X^\\U\\^ + {A,Lp) + {U,Lt). 

In the normal case A = — 1, and the maximality condition (I4.6P easily implies 
that each normal extremal curve {g{t), L(t)) is an integral curve of the Hamilton- 
ian 

(4.7) H=^\\U\f + {ALp) 

generated by the extremal control U{t) — Lt{t). This Hamiltonian will be referred 
to as the affine Hamiltonian. 

In the abnormal case the maximization relative to U results in a constraint 
if — and does not directly yield the value for U. Further investigations of these 
extremals will be deferred to the next section. 



4.1. Extremal equations. 

Hamiltonian equations (4.1) reveal their symmetries more readily when recast on 
the Lie algebras rather than on their duals. In order to treat the affine Hamiltonian 
both as a Hamiltonian on G and as a Hamiltonian on the semidirect product Gs, 
equations (4.1) need to be recast on g and Qs- Since the Lie bracket is different in 
two cases the differential equations take on different forms. 

Recall that [A, i?]^ denotes the Lie bracket that deforms the semisimple Lie 
bracket when s = 1 to the semidirect Lie bracket when s = 0. Let dh = dhp + 
dht, L — Lp + Lf,, X = Xp + Xi denote the appropriate decompositions relative to 
p and 6. 

Then §{X) = -ad*{dh{l){lit)iX) = -l{[dh,X]s) corresponds to {^,X) = 
— {L, [dh,X]s) ; the latter implies that 

{^,Xp) + {^,Xt) = -{Lp,[dhp,Xt] + [dh,Xp])~{Lt,[dh,,Xt] + s[dhp,Xp]) , 

or 

{^,Xp)+ = {[dht,Lt] + [dhp,LplXt) + {[dht,Lp] + s[dhp,LtlXp). 

Therefore, 

(4.8) ^ - [dht,Lt] + [dhp,Lpl^ = [dht,Lp] + s[dhp,Lt]. 

In the case of affine Hamiltonian H given by (|4.7p . dH = L^+A and the preceding 
equations become 

(4.9) ^(t) = [Lt,Lp]+s[A,Lt]^[sA-Lp,Lt],^it) ^[A,Lp]. 

Abnonormal extremals are the integral curves of Hq = {L, A + U) subject to the 
constraint — 0, that is, abnormal extremal curves {g{t), L{t)) are the solutions 
of 

(4.10) I = g{t)iA + Uit)), 

(4.11) ^(t) = [Uit),Lp]+s[A,Lt],^it) ^[A,Lp] + [Uit),Lt] 

subject to Lt{t) = for all t e [0, T]. They are described by the following 
proposition 
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Proposition 6. Abnormal extremal curves are the solutions of = g{t){A + U{t)) 
generated by bounded and measurable controls U{t) £ t that satisfy the constraints 

(4.12) [A,Lp] =0, [Lp,U{t)]^0 

for some element Lp in p. 

If Lp is regular, then the corresponding abnormal extremal curve whose projection 
on G is optimal is also normal. 

Proof. Suppose that g{t), Lp{t), Lt{t) = is an abnormal extremal curve generated 
by the control U{t) in {. Then equations (|4.10l) imply that = [U{t),Lp] and 
[A, Lp] = 0. This means that Lp{t) belongs to the maximal abelian subalgebra A in 
p that contains A. Therefore, also belongs to A. 

If X is an arbitrary element of A then {[U{t),Lp] ,X ) = {U{t), [Lp,X]) = 0. 
Therefore, {X, ~ 0. Since the Killing form is nondegenerate on A, — 
and therefore, Lp (t) is constant. This proves the first part of the proposition. 

To prove the second part assume that Lp is regular. Then [A, Lp] = implies 
that [Lp, [A, U{t)]] — 0. Since Lp is regular and belongs to A, [A, /7(t)] also belongs 
to A. It then follows that [A, U (<)] = by the argument identical to the one used 
in the preceding paragraph. 

It now follows that g{t) = g{0)e^^h{t) where h{t) denotes the solution of ^{t) = 
h{t)U{t), h{0) = I. Let go and gi denote the boundary points relative to which 
g{t) is optimal. Then h{t) is optimal relative to h{0) = I and h{T) = e~^^gi. This 
means that h{t) is a geodesic in K relative to the bi-invariant metric induced by 
(, ). Hence the control that generates h{t) must be constant, i.e., h{t) — e*^* for 
some element [/ in 4. 

The reader can readily verify that each trajectory {g{t), U{t)) of the affine system 
in which the control U{t) is constant and commutes with A is the projection of a 
solution of (gS)). 

□ 

Corollary 2. Ifpis such that each non-zero element is regular, then each abnormal 
extremal that projects onto an optimal trajectory is also a projection of a normal 
extremal curve. In particular, on isometry groups of space forms ( simply connected 
symmetric spaces of constant curvature ) each optimal trajectory of the affine system 
is the projection of a normal extremal curve. 

Proof. See the discussion on space forms in Section 6. □ 

Remark 1. The above proposition raises an interesting question. 

Is every optimal trajectory on an arbitrary symmetric space the projection of a 

normal extremal curve ? 

It seems that G — SLn{R) with p the space of symmetric matrices with trace zero 
and t = s0n{R) is a good testing ground for this question. In this situation there 
are plenty abnormal extremal curves but it is not clear exactly how they relate to 
optimality. 
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5. Spectral representation and its consequences 

We will now recall an observation made in ([18] ) that a system of differential 
equations of the form 

admits a spectral representation 

(5.2) ^^[A'h,L,] 

with M\ — jXq and L\ = jXq + Xi + \X2 + ■ ■ ■ X^^Xn+i- This representation 
is a consequence of a dilational symmetry 

(5.3) Xq = "J^Xq, Xi = Xi,X2 = XX2, ■ ■ ■ , Xn+i — A"X„+i. 

For then Xq, Xi, . . . , Xn+i also satisfy f lS-ll *) and therefore, ^ = [Xq, L] where 
L = Xq + Xi + ■ ■ ■ Xn+l- 

Extremal equations (|4.9p are of the form (|5.in with Xq = Lp—eA, Xi = — L{ and 
X2 = A. Since ( 15. 2p is invariant under a multiplication by A it will be convenient 
to redefine Lx as L\ = Xq + XXi + X^X2 + ■ ■ ■ X'^'^^Xn+i in which case 

(5.4) Lx=Lp- XLi + (A^ - s)A 



is a spectral matrix for equations (j4.9p , in the sense that the spectral invariants 
of Lx are constants of motion for the corresponding Hamiltonian system. Moreover, 
these functions are in involution according to the following proposition. 

Proposition 7. The spectral invariants of Lx = Lp — XLi + (A^ — 1)A Poisson 
commute with each other relative to the semisimple Lie algebra structure, while 
the spectral invariants of Lx = Lp — XLi + X^A Poisson commute relative to the 
semidirect product structure. 

The proof below is a minor adaptation of the one presented in ([Ml) and ( |2Bi ) . 

Proof. Let T : g* ^ g* be defined by r(p + fc) = - fc + for p e p*, fc € T . 
Here, a is a fixed element of p* and A and are parameters. Then, T^^ = 
Xp—k — A/ia. This diffeomorphism extends to a diffeomorphism on forms according 
to the following formula: 

(5.5) {/, g}xAi) = (T o {/, g}m = {/ ° T-\g o T-'}{TiO), 

where { , } denote the canonical Poisson form on g* (relative to the semisimple 
structure). A simple calculation shows that 

(5.6) {/, <?}a,m = -A'{/, 9} - Xf4f, g}a - (1 - X^){f, 5}, 

where {f,g}a = {/, 5}(q.) and {f,g}s is the Poisson bracket relative to the 
semidirect product structure. Relative to the semidirect structure {f,g}s the 
shifted Poisson bracket ({/, g}s)\,tj.{£,) = {To{f, g}s){0 takes on a slightly different 
form:({/,g}s)A,^ = -{f,g}s - A/i({/,g}s) a 

Functions on g* which Poisson commute with any other function on g* are 
called Casimirs., i.e., Casimirs are the elements of the center of the Poisson algebra 

C°°(g*). 
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It / is any Casimir then fx^fi — T o f satisfies {f\.^i,g}\.ti = for any 
function g on g*and any parameters A and /i. In the case that g is another 
Casimir then f\^^ and g\^^ = T o g satisfy 

for any values Ai, /ii, A2, /i2-The same applies to the semidirect Poisson bracket. 
Suppose now that = ^ . It follows from (|5.6p that {/, (7}a,^ — —^^{f,g} — 
(l-A2)({/,.g}, + {/,.9}„). Therefore, ^ 

~ A2_1 {/Ai,M1 ' 5a2,H2 i-'ll.Ml ~ A'-^ — 1 ^-^^I'/^l ' ^^2.f 2 J" -''2, A' 2 ~ (1_ A^ ) ( 1 — A^ ) •^'^1 ''^l ' ^^2 .A»2 J" ' 

Since Ai, and Ai are arbitrary {/ai.^i , gx^ ,,2^ ~ ^' This argument proves the first 
part of the proposition because AT ^ Lp — XLf + (A^ — I) A when fi — ^ after 
the identifications p — > Lp, k ^ and a ^ A. 

In the semidirect product AT = Lp — AT{ + X^A when X — ^. Then, 

= A^{/Ai,M1 I 3A2.fi2 J'") -^l.Ml ~ Af ({ Al^Pl I 5A2.fi2 J''*) >'2,A'2 " ( AJ ~ A^ ^ •^'''1 ^'^1 ' ^^^2 .f'2 J"" ' 

Therefore, {/ai.^i , 5a2,m2 J''' = ° 

6. Specific cases 

It will be convenient to single out some symmetric pairs (G, K) on which more 
detailed integrability investigations can be carried out. 

6.0.1. Non compact Riemannian symmetric spaces. Every semidirect prod- 
uct G = V ><i H admits an involutive automorphism a{v,h) = (—v,h),{v,h) G 
V t< H with K — {0} X H the group of fixed points under a. The corresponding 
splitting is given by p = ^ x {0} and t = {0} x t). The pair (G, K) is a Euclidean 
symmetric pair. Below are some examples of non-Euclidean symmetric spaces. 

Selfadjoint groups. A matrix group G is called self adjoint if the transpose g^ 
belongs to G for every g ^ G. Let G C SLn{R) be any self adjoint group. Define 
a : G G by a{g) = {g^)-\ Then a{g) ^ g if and only if g £ SOn{R) n G. 
Assuming that G 7^ SOn{R), then a is an involutive automorphism and {G,K) is 
a symmetric pair with K = GD SOn (R) ■ The splitting of g induced by a is given 

by 

(6.1) p ^ {Ae g: A^ ^ A}, and t ^ {A e g : A^ = -A}. 
The quadratic form defined by 

(6.2) {A, B) = ^Trace{AB) 

is Adfc invariant and positive definite on p and negative definite on 5. Hence, 

(G, K) is a Riemannian symmetric pair. Below are some noteworthy special cases 

of self adjoint groups. 

Positive definite matrices. G = SLn{R), K = SOn{R). Then SL,,{R) / SOn{R) can 

be identified with the space of positive definite n x n matrices with real entries. 
The generalized upper half plane. G = Spn, K = S02n H Spn = SUn- 
Recall that Spn denotes the group that leaves the symplectic form < x,y >= 

SiLi ^iUi+n — UiXi+n in IR-^" invariant. In this situation 
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The quotient Spn/SUn can be considered as the generahzed upper half plane 
since it can can be realized also as the space of complex nx n matrices Z with 

(6.3) Z = X + iY 

with X and Y real n x n symmetric matrices and Y positive definite. 

Forn = 1, Spi = SL2, SUi = S02iR) and Spi/SUi coincides with S'L2(i?)/5'02(i?). 
The latter, with its Ricmannian metric induced by ( , ) , is identified with Poincare's 
upper half plane. 

Open sets of Grassmeinnians. G = SO{p,q), K = {Op{R) x Oq{R)) fl 

SOp+q{R). 

The quotient space M+g = SO{p, q)/{Op{R) x Oq{R)) n SOp+q{R) is identified 
with the open subset of Grassmannians Grq(M.P^^') consisting of all q dimensional 
subspaces in on which the quadratic form {x, y)p,q = — XiUi + 

J2i^p_^_i XiDi is positive definite. The corresponding Lie algebra splitting is given 

by 



T 



p = {M = [ Q ] ^ P 1 fnatrix},and 

The space Mi+„ can be identified with the hyperboloid H" = {a; e M"+^ 
j — + • • • + .XjJ = 1, xo > 0} via the following identification. Let P dcnot 







the orthogonal complement relative to ( , of an n dimensional subspace Q in 
]R"+^ on which {x,y)i^n = —xoVo + X^"=2^i2/« positive definite. Since ( , 
is positive on Q and non-degenerate on P is transversal to Q and hence is 

one dimensional. Let p = (po, . . . ,Pn) be any non-zero point of P. Since the form 
( , is indefinite on M"+^, {p,p)i,n < 0. If p is normalized so that po > 
and {p ,p) = —1 then p e H" and Q is identified with the tangent space at p. 

6.0.2. Compact Riemannian symmetric spaces. 
The Grassmannians. Let G = SOp+q{R) with the automorphism (j{g) = 
JgJ"^ where J denotes a diagonal matrix with its first p diagonal entries equal 
to 1 and the remaining diagonal entries equal to —1. It follows that (j{g) = g if 
and only if J = gJg^ ■ An easy calculation shows that J = gJg^ if and only if 

^ ) where gi is apxp matrix, 52 is a q x q matrix and gi = gf and 
U 52 / 

.92 = 92 ■ Hence the isotropy group K is equal to {Op{R) x Oq{R)) fl SOp+q{R), 
which will be denoted by S{Op{R) x Oq{R)). 

The tangent map of a splits g =sOp+q{R) into p, the vector space of matrices 

P = ^ ^ ^ where B is a p x q matrix, and t the Lie algebra of K 

consisting of matrices Q ^ ^ -C* ) ^^^^ ^ ^ antisymmetric. 

The pair (G, K) is a symmetric Riemannian pair with the metric on p defined 
by the quadratic form (Pi,P2) = — 5T'r(PiP2). The homogeneous space Gvp = 
SOp+q{R) / S{Op{R) X Oq(R)) is the space of p dimensional linear subspaces in 
MP+« and it is a double cover of SOp+q{R)/ SOp{R) x SOq{R) . The latter is the 
space of oriented p dimensional linear subspaces in MP'^'^. 
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When p = 1 and q = n, then the set of oriented lines in M"+-'^ is identified with 
the sphere §" and the above gives 

= SOn+l{R)/{l} X SOn{R). 

Complex symmetric matrices. Let G = SUn with <7{g) — {g^)~^. It 
follows that a{g) = g ii and only if 5 G S'0„(J?).The corresponding splitting of su„ 
identifies 6 with the real part of matrices in sm„ and p with the imaginary matrices 
in SUn, i.e., p = {iY : Y = Y^} and t = {X : X'^ = -X}. 

The pair (SUn, SOn{R)) is a symmetric Ricmannian pair with the metric induced 
by the trace form {A. B) = -^Tr{AB). The quotient space M = SUrJSOn{R) can 
be identified with complex matrices of the form e*'^ for some symmetric real matrix 
A, because every matrix g G SUn can be written in its polar form as g = e^^R 
for some R e 50„(i?). 

For n = 2, M is a two dimensional sphere as can be verified by the following 

argument. If A = i ^ ^ ^ denote a matrix in p then = — (a^ + 6^)7 and 
therefore, 



e'^ = / cos \/a2 + 62 + A sin \/aF+¥, 



or 



_ / cos Va^ + b^ + sin V a^ + sin Va^ + b^ \ 

\^ _^sinv^^2TF cosV^^TF-^^sinV^^TF j 

Then x = cos VaF^\-W, y = ^J^f^-i sin + 6^ and z = ^^^j^yi sin Va^ + 6^, 
identifies the above matrix with the sphere + + z'^ = 1. The decomposition 
g = e^^R corresponds to the Hopf fibration ^ S"^ ^ S^. 

6.1. Space forms. Simply connected Ricmannian spaces of constant curvature, 
known as space forms, consist of hyperboloids (spaces of negative curvature), 
spheres ( spaces of negative curvature) and Euclidean spaces (spaces of zero curva- 
ture). The normalized prototypes are the unit hypcrboloid H", the unit sphere S" 
and the Euclidean space E". It follows from above that 

(6.4) S" = SO„+i/K, W = SO{n, 1)/K, E" = M" k SOn{R)/K, 

where K = {l}x SOn{R). 

The splitting of the corresponding algebras can be described in terms of the 
curvature parameter e = ±1, with 

(6.5) P^ = {(° T^)'fe^"}'«^ = {(S X )' 

It will be convenient to introduce a shorthand notation and write 

(6.6) M, = G,/K,K={l}xSOn{R), 

with Ge equal to SOn+i{R) when e = 1, SO{n,l) when e = — 1, and SEn when 
e = 0. 
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7. Affine problem on space forms 

7.1. Elastic curves and the pendulum. On space forms equations ()4.9|) admit 
additional integrals of motion in involution with the spectral ones described by 
Proposition [7l 

They are described as follows: let tA ~ {M E i : [AI,A] = 0}, and let tj^ denote 
the orthogonal complement of iA in i relative to ( , ) . It is easy to see that f ^ is a 
Lie subalgebra of 6 and that [A, Lp] £ ij^. Therefore, the projection of i{ on £^ is 
constant along the solutions of (j4.9p . 

Recall now that an extremal control U{t) is equal to Lt(t) and that the cor- 
responding extremal energy is equal to i /o \\Uit)\\^dt. Let L((t) =La + L\{t) 
denote the decomposition of L{(t) onto the factors Ia and Then, 

^ \\U{t)\\^dt= \ r {\\LAf + \\Lim)dt = \ C{\\L\{t)f+ constant 



2 1 M V/II 0/ II /IV/II/ r\ 

■10 ^ Jo Jo 



along each trajectory of (|4.9p . 

Remarkably, — K^t)"^, where K{t) denotes the geodesic curvature of 

the projected curve on G^/K, whenever La = and \\A\\ = 1 ( the norm of a 

matrix ^ = ) ^^^^"^ VJ2"=i «l )■ 

To demonstrate this fact, consider as a principal Gc bundle with connection 
T) consisting of left invariant vector fields on Ge that take values in at the group 
identity. In this setting curves g{t)m G^ are called horizontal if ^ G ^(.9(0) 
equivalently, if g~^(t)^{t) S pe for all t. Curves x{t) in G^/K can be represented 
by horizontal curves via the formula 







V / 



In this representation, ||^||e, the Riemannian length in of the tangent vector 
f , is given by \\A{t)\l where A{t) = ^{t)g-Ht). 

Solution curves g{t) of the afhne system ^(t) — g{t){A + U{t)) project onto 
the same curve x{t) as the associated horizontal curves g{t) = g{t)h{t), where h{t) 
is a solution in K of f^(t) = h{t)U{t). It follows that ^{t) = g(t){h{t)Ah-^{t)). 
Hence, 

\\'^^\\,^\\h{t)Ah-\t)\\ = \\A\\^l. 

Then, ^{^), the covariant derivative of ^ along x{t), is given by 
D fir 

-^{-)^mmm{t),A]h-\t)) eo. 

ax dt 

Since ||^||e = \\A\\ = 1, the geodesic curvature K{t) of x{t) is given by 

In particular along the extremal curves U{t) — Lt hence, = ||[i;4(t), 

when La — 0. 
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It remains to show that ||[L;4(t), = ||L;4(i)|p. There exists h E K such that 
h ^Ah^Ei— f j , a consequence of the fact that K acts transitively 

by adjoint action on the unit sphere in p^. Then, 

= h-'[A,tA]h ^ [Ei,h-^tA% 





) ^h-Hih^il -F I , ; e 

so„_i(i?) 



K^t) = \\[Li{t),AW - \\h-'L\it)h,E,\\^ = ||ei A/|p = H^lp. 
Therefore, = ||L;4(t)||, and the extremal energy i jJ'dlLj^lpdi is given by 

1 

2 X 

Curves a;(i) in the base space Gt/ K which are the projections of these extremal 
curves are called elastic and ^ K'^{t)dt is called their elastic energy f|16|). Equa- 
tions (|4.9p with = can be rephrased as 

(7.1) ^(t)^[Li,L,] + s[A,Li],^{t) ^[A,L,]. 

We will return to these equations after a brief digression to mechanics and the 
connections between the elastic problem and the motions of a mathematical pen- 
dulum. 

7.1.1. The pendulum. There is a remarkable (and somewhat mysterious) connec- 
tion between elastic curves and heavy tops that will be recalled below ( see also 
([17j) for a more general discussion). Consider first an n dimensional pendulum 
of unit length suspended at the origin of M" and acted upon by the " gravitational 
force" F = —ei, where ei, . . . , e„ denote the standard basis in M"(here, all physical 
constants are normalized to one). 

The motions of the pendulum are confined to the unit sphere S"~^. For each 
curve curve q{t) on let fi{t), . . . , fn{t) denote an orthonormal frame , called 
the moving frame, adapted to q{t) by the constraint q{t) = fi{t) and positively 
oriented relative to the absolute frame ei, . . . , e„, in the sense that the matrix R(t) 
defined by fi{t) = R{t)ei, i = 1, . . . , n, belongs to SOn{R)- 

This choice of polarization identifies the sphere as the quotient G/K with G = 
SOn{R) and K the isotropy subgroup of SOn{R) defined by Kei — e\. Evidently, 
K — \\\ Y. SOn-i{R)- Let to denote the Lie algebra of K and let denote the 
orthogonal complement in g = so„(i?) relative to the trace form. Then 

We will regard SOn (R) as the principal SOn {R) bundle ( under the right action) 
over §"~^ with a connection D consisting of the left invariant vector fields with 
values in ^i. As usual, vector fields in 2? and their integral curves will be called 
horizontal. 
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It follows that every curve q{t) on S"~ can be lifted to a horizontal curve R{t) , 

in the sense that q{t) — R{t)ei, and ^ — R{t) ^ J^^^ ^ for some curve 

u{t) in R"~^. Furthermore, it follows that any two such liftings are related by a 
left multiple by an element in K. 

The kinetic energy T associated with a path in §"~^ is given by 

l||C?gil2 l..dR 2 ^\\T:>i-\f \ \\2 1|| / M|2 



The potential energy V{q) relative to a fixed point point qq is given by V{q) ~ 
— J^^ F- ^dt, where a{t) is a path from qq to q. It follows that V{q) — ei • {q — qa). It 
is convenient to take qo = — ei in which case V = ei ■ q + 1. 

The Principle of Least Action states that each motion q{t) of the pendulum 
minimizes the action J^^ C{q{t), ^)dt over the paths from q{tQ) to q{ti) for any 
to and ti (sufficiently near each other to avoid conjugate points), where C denotes 
the Lagrangian C = T — V. Thus motions of the pendulum can be viewed as the 
solutions of the following optimal control problem on S'0„(i?) : 

Minimize the integral — (ei • Rei + l))dt over the solutions R{t) £ 

SOn{R) of ^ = R{t)Ui{t), Ui{t) = ( ^ . ^^n^*^ ) subject to the boundary 



conditions R{to) £ {R : Reo — qo}, R{ti) G {R : Reo = qi}- 

The Maximum Principle then leads to the energy Hamiltonian H on the cotan- 
gent bundle of S'0„(i?). In the realization of the cotangent bundle as the product 
SOn{R) X so* (i?), further identified with the tangent bundle SOn{R) x so„(i?) 
via the trace form {A,B) = —^Tr{AB), the energy Hamiltonian is given by 

= ^(Qi, Qi) + (ei • Rei + 1) where Qi denotes the projection on 6i of a matrix 
Q in so„ {R)- 

The Hamiltonian equations associated with H (see ([H], Ch. IV) for details) 
are given by : 

(7.2) ^ = R{t){Qi{t)), ^{t) = mt), Q{t)] - R^{t)e, A ei 

It is evident from ( 17.21 ) that the projection Qo oiQ{t) on {q is constant. But then 
this constant must be zero because of the transversality condition imposed by the 
Maximum principle. To be more explicit, recall that the transversality condition 
states that each extremal curve {R{t), Q{t)) annihilates the tangent vectors of 5*0 at 
i?(to). Therefore, {Q{to),X) = for aU X £ to { since the tangent space at R{to) 
is equal to {R{to)X : X £ to})- The transversality condition at the terminal time 
ti reaffirms that Qo — 0; hence, it and is redundant in this case. 

Equations (|7.2I ) can be lifted to the semidirect product M" x so„ (R) by iden- 
tifying vector p{t) — ~R^{t)ei with matrix P[t) = ^ ^ ^ and vector ei 

with matrix Ei — { ^ n I ■ Both matrices belong to the Cartan space p in the 

semidirect Lie algebra { f ^ X i ' ^ X£sOn{R)}. 
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Let Q ~ ^ Q Q ^ denote the embedding of so„ (R) into the semidirect product 

algebra, and let Qq and Qi denote the embeddings of Qq and Qi. 
Then 

(^•3) ( QMt) ) = t^W'^^]' "^'^ [^-^W] - ( pit)\e. 

It follows that the extremal equations ( 17.21 ) can be written also as 



(7.4) J = gitm + Ki{t)), ^ = [Pit),Qil ^ - [EuP{t)h 

where g{t) 



1 

qit) R{t) 

The Lie algebra part of equations (|7.4p agree with equations ( 17.11 ) when s — 0, 
e = 1 and A = Ei. So on the level of Lie algebras, the equations of the mathematical 
pendulum coincide with the equations for the Euclidean elastic curves. 

Consider now the isospectral matrix L\ = Lp — AL;^ + (A'^ — e)A associated 
with { 17.11 ). Since the spectral invariants of L\ are invariant under conjugations 
by elements in K there is no loss in generality if A is taken to be Ei. Matrices in 

-ep^ 



pe, e = ±1 are of the form |^ ^ y ^'^^ '^^^ written as p Ae cq where 

V = i ^ I , P € K" and where 
\P J 

n 

(7.5) {a Aeb)x ^ {b,x)^a - {a,x)eb, with {x, y)e = xqUq + ep'^ Xjifi, 

i=l 

( \ 

while matrices in can be written as x A ei with = , a; G K." ^ . It 

follows that the range of L\ is contained in the linear span of p,x,eo,ei. An easy 
calculation shows that the characteristic polynomial of the restriction of L\ to 
this vector space is given by 

(7.6) e + + C2 = 



with 



Ci=6(A2-,s) + (A2-.)H + .||Li|p + ||Lp|p,C2 = A2(||Liin|Lp|p-||[Li,Lp]|p 



Mil ^ll-^pll :'-2— --v VII-i^All II-^PII -IIL-^A'- 



The coefhcients ci and C2 reveal I2 = ||L;4|p||ip|p — ip]|p as a new integral 
of motion in addition to the Hamiltonian 7i and the Casimir Ii = s| | p + | |Lp | p. 
This integral of motion admits a nice geometric interpretation relative to the 
underlying elastic curve x{t): 

where r(i) denotes the torsion of x{t). Then ^(t) = K^(i) is a solution of 

(7.7) (^)2 = -e + m-k)e + ^{h - n^)^ + Ah 

and hence is solvable in terms of elliptic functions ([16], [17] ). 
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These geometric identifications suggest an integrating procedure in terms of the 
Serret-Frenet frames along an elastic curve Let T,N,B denote the Serret- 

Frenet triad given by the standard formulas 

(7.8) ^ = >^mit), ^ ^ -^(i)T(t) + rit)Bit), ^ = -n{t)Nit). 

The tangent vector T{t) can be identified with T{t) = {h{t)Eih-^{t)) in the 
horizontal distribution I? where h{t) is a solution of ^{t) = h{t){U{t)) with 
U{t) = L-^{t) . Then the normal and the binormal vectors N{t) and B{t) can 
be easily obtained from (|7.8|) . It was shown first in [9] and then in [16] that ^ 
is contained in the linear span of T{t), N{t), B{t). Hence, equations (|7.8p carry 
complete information about elastic curves. By analogy, the equations of the math- 
ematical pendulum are also integrable by an identical procedure. 

It can be shown that the general case with La an arbitrary constant is related 
to an n dimensional heavy top with equal principal moments of inertia, but these 
details will not be addressed here. 

Remark 2. The affine distribution T^ig) = {A + J7 : J7 G €} does not extend to 
the elastic problems on more general symmetric spaces because the isotropy group 
K does not act transitively on the spheres in the Cartan space p. Hence not every 
curve x(t) in G/K can be lifted to a horizontal curve g{t) that is a solution of 
^ = g(t){{h(t)Ah~'^ (t)), for some curve h{t) in K. This observation raises a 
question about the geometric significance of the affine problem for general symmetric 
spaces. 

8. Affine problem on coadjoint orbits 

Certain coadjoint orbits coincide with the cotangent bundles of quadric surfaces 
and the restriction of the affine Hamiltonian to these orbits coincides with the 
Hamiltonians associated with mechanical systems with quadratic potential. On 
these orbits the spectral invariants of (|5.4I ) form Lagrangian submanifolds of the 
orbits, or, stated differently, the restrictions of the Hamiltonian to such manifolds 
become completely integrable. These findings provide a natural theoretical frame- 
work for several classically known integrability results and at the same time point 
to a larger class of systems that conform to the same integration procedures. The 
text below supports these claims in complete detail. 

Recall that the coadjoint orbit Og{Io) of G through Iq e g* is defined 

OG{h) ={l ■.l = Adl-^{h).g(.G}, 

where ^rfg-i (^o)(A") — lQ{Adg-i {X)), A" e g. Also recall that g* is a Poisson mani- 
fold under the Poisson bracket {/, h}{l) — l{[df, dh]), I ^ q* , and that g* is foliated 
by coadjoint orbits of G each of which is symplectic. More precisely, the tangent 
space of Og{Io) at I consists of vectors v = ad*M{l), A/ e g and the symplectic 
form a; at Hs given by 

(8.1) uJi{vi,V2) = /([Afi,M2]), with vi = ad*Mi{l) and V2 = ad*M2{l). 
([2], Appendix 2). 

In the semisimple case Og{Iq) is identified with the adjoint orbit AdGiLo) via 
the correspondence (L , X) = 1{X) for all X g g. Consequently, each adjoint orbit 
in a semisimple Lie algebra is even dimensional. In this correspondence, tangent 
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vectors at I are identified with matrices v — [L, M] and the symplectic form takes 
on its dual form ujl{vi,V2) = {L, [Mi,M2]). 

When (G, K) is a symmetric pair than g* carries another Poisson structure { , }s 
induced by the semidirect product p x f. As in the semisimple case the quadratic 
form ( , ) can be used to identify the coadjoint orbits with certain submanifolds of 
Q. Since the Kilhng form is not invariant relative to the semidirect Lie bracket, 
these manifolds need not coincide with the adjoint orbits. The proposition below 
describes their structure. 

Proposition 8. Suppose that Iq e Q* , g — {X,h) £ Gs and I = Ad*_i{lQ), and 
further suppose that Iq — > Lq = Pq + Qo, CLnd I — > L = P + Q are the cor- 
respondences defined by the Killing form with Pq and P in p and Qq and Q in 
t. Then 

(8.2) P = AdhiPo), and Q = [AdhiPo), X] + AdhiQo)- 

Proof. Let Z = U + V he an arbitrary point of g with U E p and V E i. 
Then 

Adg-^{Z) = £{g-\eU,e^^)g%=Q = ±{~Adn-.{X)+Adh-i{eU+e^''{X)e-^-''),h-^e^^h) 

^Ad^-i{U+[X, V])+Adj,-^{V). 

Hence, 

1{Z) = h{Adg^i{Z) = (Po,A4-i(C/+ V])) + {Qo,Adh-i{V)) 
= {AdhiPo), U) + (M>(Po), [X, V]) + {AdhiQo). V) 
= {Adh{Po),U) + {[Adh{Po),XlV) + {Adh{QQ),V) = {P,U) + {Q,V) 

Therefore, 

P = Adh{Po), and Q - [Adh{PQ),X] + AdniQ^). □ 

For left invariant Hamiltonians H each coadjoint orbit is an integral manifold 
for the Hamiltonian vector field H. Moreover, the Hamiltonian vector field on a 
coadjoint orbit induced by the restriction of H coincides with the restriction of H 
to the coadjoint orbit. The latter fact will be of central importance for the rest of 
the paper as it will be shown that the Hamiltonian (|4.7p restricted to the coadjoint 
orbits through rank one matrices in SLn+i{R) relate directly to Kepler's problem, 
geodesic problem of Jacobi and the mechanical problem of Newman. The identifi- 
cation with the affine problem provides natural explanation for their integrability. 

8.1. Coadjoint orbits on the vector space of matrices of trace zero. The 

vector space Vn oi n x n matrices of trace zero admits several kinds of Lie 
algebras and each of these Lie algebras induces its own Poisson structure on the 
dual space V* . The most common Poisson structure is induced by the canonical Lie 
bracket, i.e., in which Vn as a Lie algebra is equal to sZ„(i?). The decomposition of 
sln{R) as the sum of symmetric and skew-symmetric matrices ( associated with the 
automorphism a{g) = {g^)~^ ) allows Vn to be considered also as the semidrect Lie 
algebra Symn x sOn{R), where Syrun denotes the space of symmetric nxn matrices 
of trace zero. There are also automorphisms of non-Riemannian type which induce 
semidirect products of their own. The paragraph below describes these semidirect 
products in some detail. 

Let n = p + q and let a : SLp+q{R) — >■ SLp+q{R) be defined by a{g) — 
J(((7"^)"^) J^^ where J is diagonal matrix with its first p diagonal entries equal to 
1 and the remaining diagonal entries equal to —1. It follows that <j{g) = g if and 
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only li J ~ gJg^ or equivalently, <j{g) — g and only if g G SO{p, q). The tangent 
map CT* given by cr* (A) — —jA^J induces a decomposition p(Bt where p consists 



A and D symmetric and B an arbitrary p x q matrix. The Lie algebra i is the Lie 



D antisymmetric. 

The symmetric pair {SLpj^q{R), SO{p,q)) is strictly pseudo- Riemannian be- 
cause the subgroup of the restrictions of Adx to p is not a compact subgroup of 
G/(p) ( the trace form is indefinite on p ). Nevertheless, this automorphism endows 
Vn with the semidirect Lie algebra p x t which will be of some relevance for the 
material below. 

8.1.1. Coadjoint orbits through rank one matrices. Consider first the sym- 
metric pair {SLn+i{R), SOn+i{R)) with the quadratic form {A,B) — —^Tr{AB) 
on sln+i{R)- This form is positive definite on the space of symmetric matrices p 
and negative definite on J = sOn+i{R)- Suppose that P — x x is a rank one 
symmetric matrix generated by a vector x G K"+-'^. Then Pq — x x — j^^^^I is 

in p since TrixiSSx) — J27=i(^i^ {x(^x)ei) — There are two coadjoint orbits 

through Pq, one relative to the action of Sln+i{R) and the other relative to the 
action of the semidirect product p x SOn+i{R)- 

Proposition 9. The coadjoint orbit S through Pq — xqi^xq — j—^jj-f is symplecto- 
morphic to the cotangent bundle of the projective space f"^^ in the semisimple case, 
and is symplectomorphic to the cotangent bundle of the sphere 5" in the semidirect 
case. 

Proof Let S denote the coadjoint orbit of Sln+i (R) through Pq — xo<^xo — 

If i? e SOn+i{R) then RPqR-^ = Rxq Rxo - j^^i^i^- K follows that xq can be 
replaced by ||xo||eo, because SOn+i(R) acts transitively on spheres. Therefore, 
S is dicffcomorphic to ||xo|p5'(eo (8i eo)5^^, S G S'/„+i(i?). Consider now the 
orbit through eo (8) bq. It follows that ^(eo (8 60)5^^ — Seo (g) (iS'-^)~^eo for 
any S G SLn+i{R). It is easy to verify that for any x and any y such that 
X ■ y — 1 there exists a matrix S G SLn+i{R) such that Scq — x and S'^y — eo. 
Hence, ^(eo ®eQ)S^^=x®y. 

The set of matricesja; ® y — ^^jp^I '■ x ■ y ~ 1} can be identified with the set of 
lines {(ax, ^y) : x ■ y = 1} which is symplectomorphic to the tangent bundle of 
pTi+i_ rpj^g latter is identified with the cotangent bundle via the ambient Euclidean 
inner product. 

Consider now the coadjoint orbit relative to the semidirect case. It follows from 
(18. 2p that S consists of matrices P = Adh{Pa) = h{xo) (8 h{xo) — j^i^I and 
Q = [Adh{PQ),X] — [x x,X] = Xx A x since Qq — 0. Therefore, P is a rank 
one matrix generated by x = h(xo), with h G SOn+i{R) and Q is a rank two 
antisymmetric matrix x h y with y — ^Xx. Since X is an arbitrary symmetric 
matrix of trace zero y can be any point in M""''^. The correspondence between 
(x, y) G R"+^ X M"+^ — )• X (g) X — ■'j^irf / -I- x A y is one to one provided that x-y = 0. 
Moreover, x can be any point of the sphere ||x|| = ||xo|| since SOn+i{R) acts 





22 



VELIMIR JURDJEVIC 



transitively on spheres by conjugations. Therefore S is identified with points 
(a;, y) such that ||a;|| — \ \xo\\ and x ■ y = which is the tangent (cotangent) bundle 
of the sphere S*" since the two bundles are identified via the Euclidean inner product 
inM"+i. □ 

It may be instructive to show directly that the canonical symplectic form on the 
cotangent bundle of the sphere coincides with the symplectic form of the coadjoint 
orbit. 

The tangent bundle of the cotangent bundle of the sphere is given by the vectors 
(a;, y, i, y) in R'*^"^"'^^ subject to the constraints ||a;|| — ||xo||, x-x^O, x-y~Oa,nd 
x-y+x-y — and these vectors are identified with matrices i'Six+x'SSx+xAy+xAy 
on the tangent bundle of the coadjoint orbit ( we have omitted the trace factor since 
it is irrelevant for the these calculations). The canonical symplectic form on the 
cotangent bundle of the sphere is given by 

^(x,y){{ii,yi), (±2,2)2)) = xiy2 - X2yi 

It follows from above that I G Oc{lo) is identified with L ~ x ® x + x A 
y. Then tangent vectors v = ad*M(l) at I are identified with matrices V via the 
formula {V, U) = {L, [M, U]) for all U e Qs ■ 

An easy calculation shows that V — [x x, B] + [x A y, A] + [x x, A] is the 
tangent vector at L defined hy M = A + B with A^ = -A, = B. Then, 
V = x(^x + x(^x + xAy + xAy implies that 

X = Ax, y — Ay — Bx 

It follows from (|8.1|) that the symplectic form on the coadjoint orbit is given by 
t^L{Vi,V2) = {L, [M2,Mi]). Then, 

'^x0x+xAyivi,V2) = {L, [B2,Ai] + [^2,-81] + [^2,^1]) = 

{x ^ X, [B2,Ai] + [A2,Bi]) + {xAy, [^2, ^i]) = 

A2X ■ Bix — Aix ■ B2X — {A2X ■ Aiy — Aix ■ A2y) = 

A2X ■ {Aiy - yi) - Aix ■ {A2y - 1/2) ~ {A2X ■ Aiy - Aix ■ A2y) = 

-A2X ■ yi + Aix ■ y2 = {±1 ■y2-x2- yi). 

Next consider analogous orbits defined by the pseudo Riemannian symmetric 
pair {SLn+i{R), SO{l,n)) with the Cartan space p consisting of matrices P = 
-p^ 



„ , with n an n X 1 matrix and Pn an n x n symmetric matrix. These 
P Po 

orbits will be defined through the hyperbolic inner product {x,y)-i — x ■ Jy = 
xiyi — J2^=2 ^iVi where J — diag{l, —1, ... — 1). It is easy to verify that P S p if 
and only if Tr{P) = and {Px,y)-i = (.T,Py)_i; similarly, A g so(l,n) if 
and only if {Ax,y)-i = —{x,Ay)^i for all x,y in M"+^. Thus p is the space of 
"hyperbolic symmetric" matrices. 

Definition 2. The hyperbolic rank one matrices are matrices of the form x <Si Jx 
for some vector x G R"^"'^. They will be denoted by {x ® x)-i. 

It follows that {x ® x) = (a;, u)^ix, and therefore, 

((x (8) a;)-iu, w)-i = (a;, M)_i(a;, u)_i — {u, (a; (g) x) _!?;)_ 1. 

Since the trace of (x x)_i is equal to ||x||?.i — (x,x)_i, P — (x (g) x)_i 

lko||"_i 



[n+l) 



■I is in p. 
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Now define rank two skew symmetric hyperbolic matrices Q — {x /\ y)~i = 
X ^ Jy — y Jx. An easy calculation shows that Q G so(l, n). Then 

Proposition 10. The coadjiont orbit S of the semidirect product Gg = Syin_^ k 
SO{l,n) through Pq = (xo^xq)-! --^I^jyj/ is equal to {{x(»x)^i-^-^^^^I+{xA 
y)-i : ||a;||-i — ||2;o||-i}- The latter is symplectomorphic to the cotangent bundle 
of the hyperboloid H" = {{x,y) : ||a;||_i = ||a;o||_i, = 0}. The canonical 

symplectic form uj on S is given by 

ujx.y{{xi,yi), (±2,2/2)) = (ii,2/2)-i - (±2 ■ yi)-i 
where {xi,yi) and (±2,2/2) denote tangent vectors at {x,y). 

Proof. The proof is analogous to the proof in the previous proposition and will be 
omitted. □ 

The coadjoint orbits of the above semidirect products through matrices of rank 
one can be expressed in terms of a single parameter e — ±1 with SO^ = S'0„+i(i?) 
for e = 1 and SO^ = SO{l,n) for e —1. Then will denote the vector space 

{ (^^ ^ : p e R", = Pa,Tr{Pa) = 0}. A matrix X belongs to if and only 

if it is symmetric relative to the quadratic form {x,y)e — xiyi + eJ2i=2 ^iVi- 

We will let denote the coadjoint orbit of the semidirect product = pe xi SOf^ 
through Pq = {xq ®Xq)^ — ^^^^j^ / where HxqUj = [x ,x)^. The symplectic form 
Wa;,y((ii, yi), (±2,2/2)) = (i2,2/i)e " (±1 ' 2/2)e is dual to the Poisson form 

ox oy ox oy 

9. The affine Hamiltonian system on coadjoint orbits of rank one 

Consider now the restriction oi H = ^{Li,Li) + {A,Lp) to the semidirect 
orbits Si where {A,B) = —^Tr{AB). This trace form is negative definite on the 
space of symmetric matrices and positive definite on so„+i(i?), but in the pseudo- 
Riemannian case it is indefinite on both p_i and 

In what follows it will be convenient to relax the condition that Tr{A) — 0. 
It is clear that both A and A — ^-^rr-/ define the same affine Hamiltonian since 

n+l 

{I,Lp) = 0. Then the restrictions of L( and Lp to are given by L( ^ (x A y)^ 
and Lp = {x ® x)i — -^^^f^I ■ An easy calculation shows that the restriction of 
H is given by 

(9.1) H =^\\x\U\\y\U-l{Ax,x),. 

9.1. Mechanical system of Newmann. If we replace A by —A then the restric- 
tion of the affine Hamiltonian to with e = 1 is given by H = ^(||a;|P||2/|P + 
^{x,Ax)) which coincides with the Hamiltonian of the mechanical problem on 
the sphere with a quadratic potential of C. Neumann ([21] and [22]). Then H = 
^(I|2^ll-ill2/|l-i + ^{x,Ax))-i can be considered as the hyperbolic analogue of the 
problem of Newmann. 

The equations of motion ()4.9p reduce to 

j^{x{t) A y{t)), = [-A, {x ® y),] = {Ax{t) A x{t)),, 
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and 

j^ixit)®xit)), - [{xit)Ayit)„ixit)®xit),] = \\x\\Uxit)(^yit))e + iyit)®xit)),). 
An easy calculation shows that the preceding equations are equivalent to 

(9.2) ^ = MMt), f W = -Axit) + ^iAx{t),xm - MtmMt). 

Equations ( 19.21 ) with e = 1 and ||a;|| — 1 form a point of departure for J. Moser's 
book on integrable Hamiltonian systems ([21 ). We will presently show that all the 
groundwork for integrability has already been laid out in this paper in the section 
on spectral representations. But first let us show that equations ( 19. 2p can also 
be derived in a self contained way from a "mechanical point of view" through the 
Maximum Principle of control. 

This mechanical problem will be phrased as an optimal control problem of min- 
imizing the Lagrangian ^ (||u(t)|p — {Ax,x)e)dt over the absolutely continuous 
curves x{t) on an interval [0,T] that satisfy \\x\\f: = ||a;o||e, ^{t) = u{t), and also 
satisfy fixed boundary conditions x(0) = xq and x{T) — xi. 

The constraint \\x\\f: = ||a;o||e implies that {u{t), x{t))^ = 0. The Maximum 
Principle of optimal control leads to the appropriate Hamiltonian on the cotan- 
gent bundle of the sphere ||a;||e = ||a:o||e = 1 (the hyperbohc sphere is the unit 
hyperboloid H"). 

We will use T*Se to denote this cotangent bundle. It will be identified with 
the subset of M"+^ x K"+^ subject to the constraints Gi = ||a;||e — 1 = and 
G2 — {yi^)t — 0. The Maximum Principle states that the appropriate Hamiltonian 
for this problem is obtained by maximizing 

Hq = —5(11^(^)11^ — (Ax, x)e) -I- (y, m),: relative to the controls u that are subject 
to the constraint Go = (x,w)e = 0. According to the method of Lagrange the 
maximal Hamiltonian is obtained by maximizing H=Hq -\- XqGq. It follows that 
the optimal control u occurs at u = y + Xqx. But then (w, a;)^ = (y, a;)e -I- Ao||a;||e 
implies that Aq = 0. Hence, the maximal value of Hq is given by 

(9.3) Ho = ^i\\y\f,+{Ax,xU 

The above Hamiltonian is to be taken as a Hamiltonian on T*S^; hence, r*S'e 
must an invariant manifold for Hq. Therefore, integral curves of Hq are the restric- 
tions to T* Se of the integral curves of a modified Hamiltonian 

H = Hq + XiGi + X2G2 

in which the multipliers Aiand A2 are determined by requiring that T*Se be in- 
variant for H. This requirement will be satisfied whenever the Poisson brackets 
{H,Gi} and {H,G2} vanish on T*Se. The vanishing of these Poisson brackets 
implies that Ai = -^^M = {f^t}' 

An easy calculation based on ()8.3|) yields {Hq,Gi} — ~2{y,x)e, {Hq,G2} = 
{Ax,x)^ — {y,y)e and {Gi,G2} — 2||a;||£ = 2 from which it follows that 

Ai = {Ax, x)^ - (y, 2/)e, A2 = (y, x)^. 



Hence, 



H= ^WvWi + ^{Ax,x)^ - {{Ax,x)^ - (y,y)JGi - {y,x)^G2 
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The flow of H restricted to Gi = 0, G2 = is given by 

.N dx dH dy dH , ^ , •, , ^ ^ , ^ 

Tfie preceding equations coincide witli (|9.2p . 

9.2. Integrability. Tire spectral invariants of L\ = Lp — AL( — X^A naturally 
lead to the appropriate coordinates in terms of which the above equations can be 
integrated. It will be more convenient to divide L\ by — and redefine Lx — 
— j^Lp + jLi + A. Since the trace of (x €5 x)e is a scalar multiple of the identity it 
is inessential for the calculations below and will be omitted. Then the spectrum of 
Lx is given by 

(9.5) = Det{zl - Lx) = Det{zl - A)Det{I - [zl - A)-^{-^Lp + iif)), 

A A 



with Lp — {x ® x)e and L^ — {x /\ y)^. It follows that spectral calculations can 
be reduced to 



(9.6) ^ Det{I - {zl - A)-\-—Lp + -L^)), Det{zl - A) ^ 0. 

A A 

Matrix I - {zl ~ A)-'^{-^Lp + \Li) is of the form 

(9.7) I - Rz{xi ® + X2 ® £.2) 

where xi = x, X2 = y, — ~j2Xi + jX2, £,2 = Rz — [zI — A)~^ . 

A simple argument involving a change of basis shows that the solution of equation 
(|9.6I can be reduced to Det{I — Wz) = 0, where = (wy) is a 2 x 2 matrix 
with entries wij equal to the {Rxi,£j)e (f^)- 

It follows that 




Then Det{I - Wz) = if and only if 1 - Tr{Wz) + Det{Wz) = which in turn 
implies that 1 + ^{RzX,x)t + j^{{RzX, x)e{Rzy,y)e - {RzX,y)^) = 0. 
Let 

(9.10) F = {Rzx, x), + {Rzx, x),{Rzy, y), - {RzX, y)l 

It follows from above that — Det{zl — Lx) outside of the spectrum of A if 
and only if F{z) = — A^. It is easy to verify that \im.z-y±oo zF{z) = (x,a;)e = 1 
which implies that F(z) takes both positive and negative values for any x =/= 0. 
Therefore, F is constant along any solution of (|9.2p . 

Function F is rational with poles at the eigenvalues of the matrix A. Hence, 
F{z) will be constant along the solutions of (19. 2p if and only if the residues of F 
are constant along the solutions of (|9.2p . 
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In the Euclidean case the eigenvalues of A are real and distinct since A is sym- 
metric and regular. Hence, there is no loss in generality in assuming that A is 
diagonal. Let ai, . . . , a„+i denote its diagonal entries Then, 

n 

F{z) ^ 



k—Q 



where Fq, . . . ,Fn denote the residues of F. It follows that Fk ~ liuiz^ak F{z)- 
Since F(2^ - V" -|- V" V" fV" ^MtA2 _ 



k=0 z-ak ^ l^k=0 l^j=0,j^k (z-ak)(z^aj) ^ 2^k=0 2^j=0J^k {z-ak){z-aj) ' 

Therefore, 

Proposition 11. Each residue F^ = x\ + X]j=o j>tfe ^'^^(ak-a')'' ' ^ = 0, ■ . ■ , ?t- «s 
an integral of motion for the Hamitonian system ( \9.2\) . Moreover, functions 
Fq, . . . , Fn are in involution relative to the canonical Poisson bracket in M"^'^ x 

Proof. The Poisson bracket relative to the orbit structure coincides with the canon- 
ical Poisson bracket on R"+i x ]R"+i. □ 

Remarks. Functions F^. are not functionally independent since '^^^^^ Fj. = H-'^IP = 
1. 

In the hyperbolic case the situation is slightly different because A cannot be 
diagonalized over the reals. In fact every regular matrix in Sym-i is conjugate 
/ -a 0^ \ 

to A = a , where a is a nonzero number and D is a. diagonal 
\ D J 

(n — 1) X (n — 1) matrix with distinct nonzero diagonal entries a2, . . . , 

The most convenient way to pass from the Euclidean to the hyperbolic case is 
to introduce complex coordinates 

(9.11) Vo = -^i^O +iXi),Vi = '~^i^0 -iXi),V2 = iX2,-..,Vn = ixn, 

(9.12) Wo = "^(yo + m),Wi = ~^(yO - m),UJ2 = iy2, ■.■,Wn^ iyn- 

Then 

((z - A)-^x, t/)_i = j^^{z{xoyo^xiyi) - a{xayi+xiyo)) - YTj=2 T^^iVi = 
I^^ow^o + TTTS^viWi + YJ-^^ = EjU I^^J^^J' provided that 

ao = iot and ai ~ ~ia. 

Therefore, the spectral function F{z) (I9.10p is formally the same in both the 
hyperbolic and the Euclidean case. It follows that 

n n n n 2 2 n 

(9-13) Fiz) = E ^ + E E 7 — — T - (E 



INTEGRABLE SYSTEMS 27 



The residues Fq, . . . , F„ defined by F{z) —J2k=o z-a given by 



9.14 = J2 fc = 0,...,n. 

Since is real valued for real z, Fi — Fq and each Fk,k ~ 2, . . . ,n, is real. 
It follows that 

Re{Fo), Im{Fo), Fk, k — 2, . . . ,n are integrals of motion for the hyperbolic New- 
mann problem. 

We leave it to the reader to show that 

1 " 1 

(9.15) -Fo = 00 + - ^— — 2 (("^o^fc ~ ^kVo) + iixiVk - ViXk)f{ak + ia), 

2 ^ " + "fe 

where <po = ^{x^ - xj) + i{xoXi + ^{xoyi - xiyo)'^), and 
with 0fe equal to 

^•^fc ^ a2^a^ "fc((^fc^o ^ a;o2/fe)^ - (a;^?/! - aJiyfe)^) + 2a((a;fcyo - XQyk){xkyi ~ xiyk) 
for k>2. 

9.3. Integration procedure. In the Euclidean case the integration procedure 
goes back to C.L. Jacobi in connection to the geodesic problem on an ellipsoid. 
Its modern version is presented in Moser's papers ([21], HO])- Rather than just to 
refer to the classical literature for details, it seems worthwhile to proceed with the 
main ingredients of this procedure. For simplicity of exposition we will confine our 
attention to the Euclidean sphere; the passage to the hyperbolic case requires only 
minor modifications. 

The integration is done on the manifold S defined by 

|N| = l,(x,y) = 0,F, = .x2+ (Mi^^^^cfe,fc = l,...,n + 1 

defined by the numbers ci, . . . , c„+i that satisfy — 1- 

The following auxiliary lemma will be useful for some calculations below 

Lemma 3. Let g{z) = Il^^i{z — Xk) where xi, . . . ,a;„ are any distinct n numbers. 
Then 

f{x) f{xk) fix) -A f{xk) 

—— = } — — -, and lim x—— = > — — r 

9{x) g [Xk)(x - Xk) a;->oo g[x) [xk) 

for any polynomial Junction f{z). 
Furthermore, ZLi J^) = « */deg(/) < n - I and ^Li JxA) = ^ */deg(/) 
= n — 1 and its leading coefficient is equal to one. 

Proof. Follows easily from the partial fraction expansion of D 
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Following Jacobi, the integration will be carried out in terms of elliptic coordi- 

71+1 

nates ui, . . . , u„ defined as the zeros of {RzX, x) = Yl2=i z-a^ ^'^^ each point on the 
sphere ||a;|| = 1. In addition, use will be made of the zeros of the rational function 
Ek=l 7^ ■ These zeros will be denoted by vi,. . . ,Vn and will be assumed all 
distinct. Let 

(9.17) m{z) = n^^i(^ - u,), a{z) = n^lHz - a^), b{z) = Ul^,{z - v^). 
Lemma 4. 

y^^ = ^, and y^^^bjz), 
f^^z-Uk a[z) t^i^~'^k a{z) 



Proof. The fact that any rational function is determined up to a constant factor 

z-cik a{z) 



by its zeros and poles implies that X^fc^i z-a ~ "^■qT?^ ' '^here c is a constant. It 



follows from Lemma [3] that c = 1 , because 

fe=o fc=o ^ ' 

The same argument carries over to 'y ?^} — □ 

Since xi = , Xk can be recovered up to a sign from ui, . . . , u„. 

Recall now the function F{z) = (1 + {Rzy,y)){RzX,x) — {RzX,y))'^. Then 
F{uk) = -{Ru.x, yf = ^^^0 ^ = ^ and therefore, 



{Ru^x.y) = ±y-^|^, fc = 1, . . . ,n. 

Each choice of the sign defines a set of n linear equations for the variables ?/i , . . . , yn+i, 
which together with (a;, y) = determine y uniquely in terms of u = (ui, . . . , u„), 
and each choice identifies . . . , m„) as a system of coordinates for the Lagrangian 
manifold S. 

Proposition 12. Vectors . . . , q^^^^ form an orthogonal frame on the 

sphere \\x\\ = 1 , with\\-§^^\\^ = "i^fe^- Suppose that {R^^x.y) = k = 

1, . . . ,ri. Then, the differential equation ^ — y is given on S by 

^TO'(ufc)^^ = ^/~a{uk)b{uk),k = 1, . . . ,n. 
The preceding equations can also be written as 

(9-18) E ,, 7 / ir-'^-^-'----'- 

t^i 2V-a(ufe)6(ufe), at 

Proof. An easy logarithmic differentiation of xl — ^ yields §^ = — — r. 
This implies that 

^^2; 1. , 

Therefore, (|^, - - K/ - A)-!^:) - 

- A)-\u,I - A)-'x, x) = - - {ukl - A)-^x, x)) = 
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4iu,-u,) ii^^j^^^) ~ i^u,x,x)) = for k^j. 
For j = k, 

-i^^^. Since T^O'l^'li;'---'^ fo™ an orthogonal frame on S. 

Let Pi , . . . , P„ denote the coordinates of y relative to the frame , , . . . , 
for G 5. It follows that = Pfc|||f^|p. But = i(K/ - 

A)-ix,y) = i(i?„,x,y) = -iyqgr 

Therefore, 

2 



m'(uk) 

Suppose now that {x{t),y{t)) is a curve in 5 with |f = y. Then, ^ = ^^^^ ^ 
ELi^fcl^- Therefore, im'K)^^ = ^ -a{uk)b{uk). 

Second expression follows from Lemma [3] which implies that X]fc=i rrJ{u ) ^ 
(^n-i.m, m < n — 1 by taking /(z) = and g{z) — m{z). Then, m!{uk) — 
2yJ —a(uk)b{uk) / which, after the substitution, leads to 

^'==1 2V-a(«fc)6("fc), ' 

J. Moser points out that equation (|9.18p is related to the Jacobi map of the 
Riemann surface 

(9.19) iv^ = -4a(z)&(z). 

In fact he shows that the Jacobi map given by 

tr[J{o,o) 2^-a{z)b{z) 

takes the divisor class defined by (u^, 2-\/a(ufe)6(Mfc)), k = l,...,n into a point 
s G C^'/T where T denotes the period lattice of the differentials of the first kind 

m)- 

10. Connection to geodesic problems on quadric surfaces: Knorrer's 

transformation 

Jacobi's geodesic problem on an ellipsoid § = {a: e M"+^ : {x,A^^x) — 1} 
consists of finding curves in § of minimal length, relative to the metric inherited 
from the ambient Euclidean metric in R""*"^, that connect a given pair of points 
in §. Jacobi was able to show that the curves of minimal length can be obtained 
from the solutions of a first order partial differential equation, known today as 
the Hamilton-Jacobi equation; in the process, he discovered an ingenious choice of 
coordinates on the ellipsoid, known today as elliptic coordinates, in terms of which 
the associated partial differential equation becomes separable with its solutions 
given by hyperelliptic functions. 

Alternatively, the geodesic equations can be represented by a Hamiltonian system 

(10 1) —-V ^ _ _ bl^lV) 4-1^ 

dt dt^ lu-ixip"^ ^- 
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on the cotangent bundle of § realized as the subset of E""*"^ x M"+^ subject to 
Gi ~ {x, A^^x) — 1 = 0,G2 — {p, A^^x) = 0. Moser shows that the above equations 
are generated by the Hamiltonian 

in R"+i X M"+i in the sense that, 

dx _dH dp _ dH 
dt dp ' dt dx 

but constrained to Gi = G2 = and H = ^, i.e., to ||p|| = 1 f|21|). 

Let us modify above equations by replacing the Euclidean inner product (z, y) by 
the inner product {x,y)e that encompasses both the Euclidean and the hyperbolic 
inner product. Then, equations (jlO.ip take on the following form: 

{10-2) —r=P, —r— — Trz — ^ — tt^A x. 

^ ' dt ^' dt \\A'^x\\^ 

assuming that ||v4~ix|| ^ 0. It follows that Gi = (x,^"^^)^ - 1 = 0, G2 = 
{p, A~-^x)e = 0, ||p||e = 1 is an invariant set for ()10.2p . 

Remarkably, equations (|10.2p can be transformed into the equations of ()9.4p by 
a transformation discovered by H. Knorrer in ([12j), and the integrals of motion 
of the geodesic problem can be deduced from the integrals of motion associated 
with the mechanical problem on the sphere. In what follows we will consider the 
inverse of the Knorrer's transformation and show that the integrals of motion for 
the geodesic problem can be deduced from the mechanical problem of Newmann. 
For that reason we will begin with equations (j9.4p written as 

nil nil 

(10.3) T=v, ^ = ~A-^u + {{A-^u,u),-\\vWi)u, \\u\U = l. e = ±1. 
as ds 

It is important to keep in mind that the matrix A also depends on e since it 
belongs to pc (modulo the trace). In what follows it will be necessary to assume 
that {Au,u)^ > 0,u^0. 

Recall that Fq — {{v,Av)f. — l){u,Au)f^ — {u,Av)'^ is an integral of motion for 
(jl0.3p as can be easily seen from (|9.10l) with z = 0. Let ^(A, u, v) — {x, p) denote the 
mapping from the manifold Nq = {(A, u,w) : ||M||e = l,{u,v)e = 0,Fo = 0, A € R} 
given by 

(10.4) a;^ p= M Av^ ) ' [ Au). 

It follows that (z, yl"-'^a;)e = 1 and that [p^A^^x)^ = 0, hence $ maps into 
the tangent bundle of the quadric {x^A~^x)^ = 1 . We will show now that $ is 
invertible and that its inverse is given by: 



10-5 A = II — , u = II II , V = (A p-{u,A p),u). 

\\A ^x\U \\A ^x\U y/{A-~^p,p), 

It follows from mM that i^(A~^p,p), ^ (^Av~\^^Au) [v- \^^u) = 
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The constraint ((v, Aw)c(it, — Au)^ — [u^Au),^ implies that {^A~^p^'p)^ 
1. Further constraints \\u\\^^ ~ 1 and = imply that {Au^u)e = ||^_\^||2 

Let m(s) and d(s) be any solutions of (110.3^ and let A(s) be a solution of 
(10 6) ^ ^ 2i^Mfll^^MliA(s) 

C?S (Au(s), u(s))e 

Then, 

^ = A^^Au = ii;, and 

^ = ^-A )(At;-|4H^A^) + ^ ±{Av-^T^^Au) = 

^{Au.u), ^ (Au,u), I ^(Aum), ds\ iAu,u), I 



^(Au,«), yJ{Au,u), ^ (Au,u)l ' ^{Au,u)^ 

Hence, 

^''■'^ Ts = r'd-s = -'^ ^- 

Now identify with a function ^. It follows from (|10.5I) that = = 

(li-ip^p)' • Hence, ||A-ix^p = '^(*) ^^'^ hence equations (|10.7p coincide with 
equations (|10.2p after the reparametrization t = 0(s) with ^ = 

The integrals of motion obtained for the mechanical problem of Newmann have 
their analogues for the problem of Jacobi via the following proposition 

Proposition 13. Let F{w) = (1 + {RiuV,v)f){RwU,u)e — {RwU,v)1, with = 
{wI-A)-'^. Then, 

F{-) = II mo! + {SzX,x))^){SzP,p)e - {SzX,p)l,Sz (z - 

z \\A ^x\\i{A ^p,p)e 
under the substitutions given by formulas \10.5\) . 

Proof. It is easy to verify that F{w) is invariant under the change of variable 
V ^ V + au with a an arbitrary scalar. Hence, F{w) = (1 + {R^V, V)){RwU, u) — 
{RwU, VY , where V = In the proof below we will use the identity 

{- - A-^)-^ +A = -{z- A)-^A^. 
z 

Then, 

1 + {RiV, V),^l + {\- A)-^)-W, V), = 

1 ((7 4\-i 42 ..^~\P.. \ ( p. ..^~\P.. \ — 

1-{{Z- A) A , A||A-lx||Je - {x\\A-ix\\,^ \\\A-^x\\J^ " 
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Further, 
and 

{Riu, V)e = -((z - A) ^A"^ \\A-ix\U^ xiti-^lwJ'^ ~ '^\\A-^x\\, ' A||l-i^|u)« = 

A||A-i£c||2 {SyX,v). 

Hence, 

^(|) = A2||A-l. || f ((l + {S,X,x)),{S,V,V)e - {S.X,p)l) ^ \\A-^x\\il{A-^p.,p)M + 
{SzX, x)e){SzP,p)e - {SzX,p)l). □ 

Corollary 3. Function G{z) = (1 + {SzX,x),){SzP,p)t — {SzX,p)1 is constant 
along the geodesic flow 110.2\) . 

Proof. \\A~^x\\'^{A~^p,p), is an integral of motion for (jl0.2p because 
i\\A-^x\\UA-^p,p),^ 

2{A-'p,A-'xUA-'p,P)e - \\A-'x\\UA-'p, (A-V,p).^n^)e = 0. 

It follows that G{z) is constant along the solutions of (|10.2p . since F{j) is constant 
along the solutions of (|9.2p . □ 



Remark 4. Function \\A ^x\\'j{A ^p,p), is known as Joachimsthal's integral of 

motion( [g;, (^). 

In the Euclidean case the matrix A can be assumed diagonal with ai, . . . , a„+i its 
eigenvalues. An argument identical to the one used above shows that 

n+l „ 

and that the residues are given by 

ri+l / \2 

(10.8) G,^pI+ ^T~"'? ,fc = l,...,r. + l, 

as reported in (|21j. The hyperbolic case differs only in minor details due to different 
canonical structures of A. 

11. The case A~0 and the problem of Kepler 

Consider now the Hamiltonian H = ^{L(,Le) on coadjoint orbits of pe xi t^, 
e = ±1, through rank one matrices in p,. We will continue with the notations 
of the last two sections and consider the coadjoint orbis through matrices Pq = 

{xq '^xo)e — ^^n+i' ^' £ = ±1 • We have seen that the these coadjoint orbits consist 
of matrices 

ip = (x (g) a;)e,i{ = (a; A y)e, = ||a;o||o {x,y), = 0, 

that are symplectomorphic to the cotangent bundle of the "sphere" {{x,y) : \\x\\e — 
llxolle, {x,y)t — 0}. The Hamiltonian equations (|4.9p and (|9.2p reduce to 

(11,1) i?|.|i.,i,|,^il.o, 
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and 



dx 1 1 1 1 2 '^y 



(11.2) 

It follows that the solutions satisfy 

'^ + \\A\l\\y\\lx^Q 

The restriction of the Hamiltonian H to these orbits is given by H = i ||a;||e||2/||e- 
Hence, on energy level H = ^, the preceding equations reduce to 

(11.3) + 

It follows that the solutions of (|11.3p are given by great circles x{t) = acos{t) + 
bs'mt for e — 1, and great hyperbolas x{t) = a cosh t + bsinht for e = — 1, with 
M', + \\b\\l = \\xo\\l (a,6)£ = 0. 

Recall now the Hamiltonian E — — associated with the problem 

of Kepler in the phase space {(g,p) S R" x R" : q ^ 0} corresponding to the 
normalized constants m = kM = 1 and the associated equations of motion 

dq _ dp 1 

dt ^' dt \\q\ 

Below is a summary of the classical theory connected with Kepler's problem. 

(1) L = qAp and F ~ Lp~ are constants of motion for (|11.4p . L is an n 
dimensional generalization of the angular momentum p x q. Its constancy 
implies that each solution remains in the plane spanned by q{0) and p(0). 

The second vector is call the Runge-Lenz vector or, sometimes, the ec- 
centricity vector. It lies in the plane spanned by p and q. 

(2) Let =2||L||2i; + l, where = -^Tr{L^) = \\q\\^\\p\\^ - [q ■ p)^ . 
Then, \\F\\ < 1 whenever E < 0, \\F\\ — 1 whenever E — 0, and ||i^|| > 
1 whenever E > 0. 

(3) A solution {q{t),p{t)) evolves on a line through the origin if and only if 
q{0) and p{Q) are colinear, that is, whenever L = 0. In the case that L ^ 

urn 



(11-4) ^^=pr^ = -jr^,. 



1 + \\F\\ COS (j){t)' 

where (j>{t) denotes the angle between F and q{t). Therefore, 

q{t) traces an ellipse when ||_F|| < 1, a parabola when ||i^|| = 1 and a 
hyperbola when ||F|| > 1. 

There is a remarkable connection between the solutions of Kepler's problem and 
the geodesic flows on space forms that was first reported by V.A. Fock in 1935 
in connection with the theory of hydrogen atom ( [8] ) which then was rediscovered 
independently by J. Moser in ([1^) for the geodesies on a sphere. Moser's study 
was later completed to all space forms by Y. Osipov in f|23)). 

As brilliant as these contributions were, they, nevertheless, did not attempt any 
explanations in regard to this enigmatic connection between planetary motions and 
geodesies on space forms. This issue later inspired V. Guillemin and S. Sternberg to 
take up the problem of Kepler in ([lOj) in a larger geometric context with Moser's 
observation at the heart of the matter. 
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It seems altogether natural to include Kepler's problem in this study. In this 
setting Kepler's system is recognized within a large class of integrable systems and 
secondly, the focus on coadjoint representations provides natural explanations for 
its connections to the geodesic problems. Following Moser we will consider the 
stereographic projection from the sphere ||x||^ = into M" given by 

Afa; — heo) + heo = (0,p), where A = . 

h — xq 

Here, {xo,xi, . . . , x„) denote the coordinates of a point x in R"+^ corresponding to 
the standard basis eo, . . . , e„. It follows that 

Consider now the extension of this mapping to the cotangent bundle of | |a;| |^ = h'^ 
that pulls back the canonical symplectic form in R" x M" onto the symplectic form 
of the cotangent bundle of ||a;||g = h^- It suffices to find a mapping q — 'i>{x,y) 
such that 

n n 

(11.6) ^ Qidpi = (y, dx)e = yodxa + £ ^ Vidxi 
1=1 1=1 

because the symplectic forms are the exterior derivatives of the preceding forms. 
It turns out that such an extension is unique by the following arguments. 
Let X = $(p) denote the mapping given by (|11.5I V Then, 

, /^"^N J / 4£/i3 2eh^ , Aeh^p-dp . 

(11.7) dx = {-Up ^ ( ^ . dp, ^ ^^^ dp - ^ 

It follows by an easy calculation that 

(11.8) IMxIl^ = dxl + eJ2dx^ = ^^^-i^llrfplp. 

Since {dx,dx), = ((O^)^dp, {§Up), = ((||)*(||)^dp, dp), = ^^^-^l^\\dp\\^ 

it follows that i^)t^ = (||p||2^gfe2-)2 In, where /„ denotes the n dimensional 
identity and (^)e the adjoint operator of ^ relative to the inner product 
( , )£,i.e.,(||)* = i^VJe, with (ll)"^ the transpose of || and diagonal 
matrix with its diagonal entries (!,£,£,. ..,e). 
Then, 

q-dp^ {y- dx)e = {y, ^dp)^ = {^)ly, dp)^ implies that q = {^)ly or, 

— w — (5^)^'^' 

since (f):y = £nMf^(||)*(||),^,. 
Equations (|11.7p reveal that 

(11.9) y=i-q.p,^^q-—p) 
from which it follows that 
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To pass to the problem of Kepler, write the Hamiltonian H = ■^\\x\\'^\\y\\'^ in 
the variables {p,q). It follows that H = Ife^e 'llPH^+f = ig (INI^+f | |g| 
The corresponding flow is given by 

(UU) ml+£fl dq__dH_ J\p\\' + eh' , 

^ ' ds dq Ah^ ^'ds dp " 2h^ ^ 

On energy level H = ijp-, ||g|p = I and the preceding equations 

reduce to 

^''■'^^ d-s='hm^^ -ds=-'i^p- 

The preceding equations coincide with the equations of Kepler's problem (|11.4p after 
the reparametrization by a parameter t = —j^ Jq \\Q{T\\dT. For then, ^ = "'R?!! 
and equations become 

dp dp ds q dq dq ds 

dt ds dt ||g|P'ds ds dt ^ 

Since ^Mte!)!||,||2 = l, 

^ - ^ll^ll' - M ^ ^(IWni.ll - 2) . ^(2 - eh^Ml - 2) . -\ek^. 

So E < in the spherical case and £' > in the hyperbolic case. 

The Euclidean case E = can be obtained by a limiting argument in which e is 
regarded as a continuous parameter which tends to zero. To explain in more detail, 
let x{t) = x{t) — xo{t)eo where x{t) is a solution of ()11.3p . If 

w{t) = lim^^o 7;Tj(2;(i)), then w{t) is a solution of 

that is, w{t) is a geodesic corresponding to the standard Euclidean metric. It then 
follows from (|11.5p that lim^^o xo ~ h and 



Moreover, lim£_j.o dxo ~ and hme^.o — dw — y^dp — j^^yP as can be 
seen from pi. 71 ). Therefore, = jj^^lMplP 



The transformation p ^ w with w — nAnP is the inversion about the circle 
IIpIP = 2 , and = jj^JiprlMplP is the corresponding transformation of the 

Euclidean metric ||rfp|P- The Hamiltonian Hq associated with this metric is equal 

toiMr||,||2. 

This Hamiltonian can be also obtained as the limit of (■^)^ ^^^^^ 4^:f^' | |g| P when 
e — ^ 0. On energy level iJ = i, ||p|p||q|| = 2 and therefore, E = 0. 

The integrals of motion for the problem of Kepler are synonymous with the 
constancy of the matrix {x A y)e along the flow of To be more specific, let 

X = xqeq + X and y = j/oeo + y- Then, 

{x A y)e = {xoeo + xA yoco + y)e = a;o(eo A y)e - yo(eo A x)^ + (x A y)e- Since 
{x{t) A y {t))e is constant along the flows of (jll.ip both a;o(eo A y)e — 2/o(eo A x)e 
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and {x /\ y)e are also constant. But then the angular momentum L = q A p and 
the Runge-Lenz vector F = Lp — are given by 

(11.13) L ^ {y A x)e and F = h{yo{eo A x)e - xqi^co A y)e)eo- 

The first equality is evident fi-om equations ( 111.51) and (|11.9|) and the fact 
that [y A x)^ = e{y A x). Second equality follows by the calculation below: 
(yo(eo A x)e - XQ(eQ A y)e)eo = -yox + xoy = 

~\h1 P)\\p\\i+sh'^P+\ MlHsh'i )( — — l-T^P) - h{-{<l-P)P^ 2 1l- 

Since ^^^t£^||g|| = 1, e/j' - ^ - IblP, and therefore, 

h{yo{eoAx)e - a;o(eo A y)e)eo = {- {q ■ p)p + WpW^ q ~ = {q A p)p - ^ F. 

11.0.1. Conic sections and the geodesies. The geodesies of the spaces of constant 
curvature are transformed into the conic sections of the problem of Kepler, a fact 
well known in conformal geometry. For the convenience of the reader not familiar 
with these facts and also for the completeness of the presentation we include the 
basic details. 

In the spherical case, the great circle x = acoswt + bsinujt with ||a|| = \\b\\ = 
h, a - b = can be rotated around eo so that a and b are in the subspace spanned 
by Co, ei, 62- Moreover, such a rotation R can be chosen so that Ra — hei. 

Let a denote the angle that the great circle makes with the plane xq = O.Then, 

Xq — h sina sin(/it), xi = hcos{ht), X2 — — /icos a sin(/i<), x^^O, i = 3;...,n + l, 

because — h? . Furthermore, {x — heo) + heo = (0,p) implies 

that 

P^ ^ ■ ^ ■ /■ ,N cos{ht), p2 = : — ^. cos a sm{ht), = 0, i = 3, . . . , n, 

1 — sm a siTi(nt) 1 — sm a sm[rit) 

and y{t) = 7^^ implies that 

?/ = (sinQ;cos(/it), — sin(/ii), — cosQ!COs(/ii), 0,..., 0). 

Then, 

IIpIP + h'^ — : . and — — — sin a cos(/ii), 

1 — sma sm[htj h 

(implied by equations (jlLOp ). Hence, 

^ 31 sinacos(/it) ^QgC^^^ 

1 — sin Q sin(/it) 1 — sin a sin(/if ) ^ ' 

V2 = 1 — . ^ . + — : ^Sttt cos a smint). 

1 — sin CK sin(/ir) 1 — sin a sin(/ir} ^ ' 

It follows that — sin(ft,i) (1 — sin a sin(ft,i)) = qi — smacos^ht) cos{ht) and 
— cos a cos{ht) (1 — sin a sm{ht)) = 92 + sin a cos{ht) cos a sm{ht) 
and therefore, 

qi = — sin(/i<) +sinQ;, q2 = — cos a cos(/it). 
Hence, the great circle is transformed into the ellipse 

{qi - sina)2 + — ^ ^ 
cos^ a 

This ellipse degenerates into a line through the origin when a = ^ , or when the 
great circle passes through heo- 

A similar argument shows that the hyperboloid x{t) = asmh{ht) + bcosh{ht) is 
transformed into the hyperbola — (gi — sina)^ H h—q2 = 1- 



INTEGRABLE SYSTEMS 



37 



In the Euclidean case, the hne w — a+bt is transformed into the curve ■fP(f)jp'P(0 
via the mapping w — j^^P- Hence, 

b-W = MF^- m^^P' Aft^'- substitutions, (p, = (u;, 

and ^ = -^^^^9, = IMF' above equation becomes 

\\'w\Y at 



Hence, 



q^]:{h{\\a\\^-\\h\\H^)-2a{{a,h) + \\h\\H)). 



2 

This equation is a parabola in the a, h plane. 

12. Concluding remarks 

The above exposition could be viewed as a first step in unifying various frag- 
mented results in the theory of integrable systems. The fact that much of this 
theory is related to Lie groups and the associated Lie algebras has been recog- 
nized in one form or another for some time now ( [3] |24j |26j |27| ). However, 
in contrast to the cited publications, the present study uses control theory and 
its Maximum Principle as a point of departure for geometric problems with non- 
holonomic constraints which greatly facilitates passage to the appropriate Hamil- 
tonians and which, at the same time, clarifies the role of the Hamiltonians for the 
original problems. 

Additionaly, the ubiquitous presence of the afRne problem on any symmetric 
space paves a way for new classes of integrable systems, for it seems very likely 
that the afRne Hamiltonian is integrable on any coadjoint orbit ( [3]). Further 
clarifications of this situation would be welcome additions to the theory of inte- 
grable systems. Along more specific lines, the study of Fedorov and Jovanovic 
(([Z])) strongly suggests that the problem of Newmann on Steifel manifolds can 
be seen also as the resticition of the affine Hamiltonian to the coadjoint orbit of the 
semidirect product through an arbitrary symmetric matrix. It would be instructive 
to investigate this situation in some detail. 

It might be also worthwhile to mention that the solutions of the affine problem on 
the unitary group would find direct applications in the emerging field of quantum 
control ( 4 ). This topic, however, because of its own intricacies is deferred to a 
separate study. 
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